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We establish in a combination of ab initio theory and experiments that the tunneling process 
in scanning tunneling microscopy /spectroscopy on the A-122 iron pnictide superconductors - in 
this case BaFe2-xCo cc As2 - involve a strong adatom filtering of the differential conductance from 
the near-^i? Fe3d states, which in turn originates from the top-most sub-surface Fe layer of the 
crystal. The calculations show that the dominance of surface Ba-related tunneling pathways leaves 
fingerprints found in the experimental differential conductance data, including large particle-hole 
asymmetry and an energy-dependent contrast inversion. 
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The discovery of high-temperature superconductivity 
in the iron pnictides has prompted a tremendous effort di- 
rected towards understanding the nature of superconduc- 
tivity in these systems. Theory indicates unconventional 
superconductivity but with a pairing symmetry that can 
be very sensitive to electronic and structural details of 
the system [1]. Experimentally, surface sensitive probes 
of the electronic structure such as angle-resolved photoe- 
mission spectroscopy (ARPES) [2-5] and scanning tun- 
neling microscopy /spectroscopy (STM/STS) [6-12], play 
a major role in the search for a microscopic understand- 
ing of the superconductivity in these systems. In this 
context, the alkaline-earth-metal AFe2As2 (A-122, A = 
Ca, Ba, Sr) iron pnicties have been widely studied due 
to the availability of high-quality crystals, and the pos- 
sibility to both n- and p-type dope these systems. It 
is therefore essential to definitively establish the proper- 
ties of A-122 surfaces. Accessing the underlying physics 
representative of the bulk relies on a quantitative under- 
standing of surface-related matrix elements coupled to 
potential surface-related changes in electronic/structural 
properties and this is exactly what we aim to do here. 

Perhaps the most prominent theater for the sur- 
face/bulk debate is the case of BaFe2-xCo x As2 (Ba- 
122) [6]. Due to the strong covalent bonding within the 
AsFe2As block, the surface of cleaved Ba-122 is expected 
to be either an exposed FeAs layer or a half layer of Ba 
atoms. Although there are indications for the former 
[12], there is a much larger body of evidence for the lat- 
ter [5, 7-10]. Electrostatics would dictate a half Ba layer 
as the crystal termination. Both the strong temperature 
dependence of the reconstructed surface and the large 
variety of topographic images observed by STM [6, 7] 
are most easily understood in terms of relatively loosely 
bound Ba atoms (which could be considered adatoms) 
constituting the termination layer. Density functional 



theory (DFT) data point to a Ba termination with a 
V2x\^2 or 2x1 reconstruction as energetically favorable 
[13]. Recently, a combined LEED / ARPES /DFT study 
showed that the experimental LEED I/V curves could 
only result from a reconstructed Ba surface and further- 
more that the associated surface-related bands overlap- 
ping with the bulk states have a strong impact on the 
breadth of spectral features observed in ARPES for this 
system [5]. 

A quantitative understanding of tunneling processes in 
complex, layered crystals is important, as is illustrated in 
the cuprate high T c systems, where asymmetry between 
electron removal and addition in STS data is seen as a 
sign of strong correlation effects [14], as well being at- 
tributed to the effect of the specific tunneling pathway 
from the Cu02 plane LDOS through the SrO and BiO 
layers to the STM tip [15]. Given the marked 2x1 and 
V2x\^2 surface reconstructions observed in the LEED 
and STM experiments on the A-122 pnictide family [6], 
it is important to uncover the impact of the real sur- 
face structure on spectroscopic tunneling measurements 
themselves. 

We tackle this issue here by theoretically establishing 
the precise role of the Ba surface in STS of Co-doped Ba- 
122 superconductors, showing that the surface Ba atoms 
strongly filter the tunneling current, leading to marked 
particle-hole asymmetry in the tunneling conductance 
dl/dV and bias-dependent contrast inversion in dl/dV 
maps. A direct comparison with experiment explicitly 
confirms these key predictions. 

We start by examining the tunneling conductance in a 
simplified Fe band-based model, incorporating the five- 
band tight-binding model of Graser et al. [16] and com- 
paring the results to the experimental data. This model 
was obtained by fitting to DFT calculations and so an ad- 
ditional two-fold renormalization of the bands is included 
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FIG. 1: (Color online) Comparison of calculated and mea- 
sured dl/dV spectra in Co-doped Ba-122. The model spec- 
trum (black) was obtained from a five-band tight binding 
model for Ba-122 including the effects of interband impurity 
scattering within the superconductor. The experimentally 
measured spectrum (red/dashed line, Vuas — —300 meV, 
I set = 50 pA) was recorded at 4 K on an overdoped sam- 
ple (x = 0.22, T c = 12 K). 



to account for electronic correlations. As a next step, su- 
perconductivity is introduced at the BCS level, including 
interband impurity scattering. Since the interest here is 
in dl/dV over energy scales almost 100 times that of A sc , 
the precision of the theoretical description of the super- 
conductivity is non-critical [17]. Fig. 1 compares the 
model dl/dV to measured spectra obtained on an over- 
doped sample (x = 0.22 i.e. 11% doping, T c = 12 K). 
On the occupied side - negative sample bias - the gross 
structure is reproduced by the model, and one might as- 
sign the broad peak in the STS trace at ~ 200 meV to 
features in the Fe-based DOS. However, the unoccupied 
states raise doubts as to whether such a simple interpre- 
tation is tenable, as the measured dl/dV exhibits a large 
particle-hole asymmetry with a much larger signal at pos- 
itive bias voltage, a characteristic wholly lacking in the 
data of the tight-binding model. This discrepancy can- 
not be removed by including additional self-energy effects 
as they would only serve to further broaden the spectral 
features. It is clear that the theory needs to leave the 
confines of the five Fe-band basis, and the remainder of 
this letter is devoted precisely to this, showing that rec- 
ognizing the filtering effect of the surface Ba states on 
the tunnel current removes the discrepancy so visible in 
Fig. 1. 

To model the tunneling process correctly we have per- 
formed DFT calculations to simulate tunneling into var- 
ious finite slab geometries of 8.5% (x = 0.17) Co-doped 
Ba-122. These calculations use the full potential local 
orbital code (FPLO), [18], in the scalar-relativistic and 
local density approximations (LDA) [19]. The choice of 
functional is not critical, since we are considering the 
non-magnetic state and fixed atomic geometries. See [17] 
for further computational details. 

The Ba-122 surface is modeled using a slab struc- 



ture with a 2x1 reconstructed, half Ba termination layer 
with the atomic positions fixed to LEED data [5]. Since 
STM/STS experiments are often performed with Pt/Ir 
tips, we model the tip with a single Pt atom (hence- 
forth referred to as the 'tip'), positioned at different lo- 
cations above the sample. The tunneling conductance 
is calculated from a weighted overlap between the slab 
wave functions and the localized 6s and 5d states of the 
tip atom: 
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where \x\) is the discrete tip state i with energy e-, and 
H is the Fermi level of the slab. The 5 function repre- 
senting the tip density of states (DOS) is replaced with 
a Lorentzian broadened expression, mimicking band for- 
mation in a real macroscopic tip. Here we used 7 = 1 
eV for all tip states, noting that the exact 7 value has 
no qualitative influence on the results, and merely results 
in a scaling of the magnitude of the resulting ^ curves. 
Eq. (1) is essentially a version of the Todorov-Pendry 
approximation [20] as has been used in a wave function 
context [21] and a Green's function context [22]. It goes 
beyond the commonly used Tersoff-Hamann model [23], 
which is derived for s-wave tips only. 

The matrix element in Eq. (1) in its full form is 
M(k,n) = (*k,n|Vr I X*), with V T the potential at the 
tip. However, our tests have shown that for a nearly 
spherical potential the effect of including V T also only 
results in a global rescale of the curves, and thus the 
form shown in Eq. (1) is sufficient. 

We first consider the matrix element for tunneling be- 
tween the tip and the sample M(k, n) = (^k,n\Xr)- Fig. 
2(a) plots the band structure for the 2x1 reconstructed 
slab, with numerous surface bands from reconstruction- 
induced back-folding and atomic distortions in the out- 
ermost As-Fe2-As block [5]. The thickness of the or- 
ange (light) shading codes for the weight associated with 
the surface-Ba character of the bands. Most of the Ba- 
dominated bands lie above E Fl except for one band 
with Ba-5d character dipping below Ep between Y and 
X. The magnitude of the tip-sample matrix element 
|M(k, n)| 2 , coded by the weight of the blue (dark) shad- 
ing, is also shown in Fig. 2(a). Comparing the two 
shadings, it is clear that the largest tunneling matrix 
elements coincide with the presence of bands with strong 
surface Ba character. To illustrate this important find- 
ing, panels (b)-(d) of Fig. 2 show the Fermi surface of 
the slab. In panel (b) the blue shading codes the sub- 
surface Fe 3d character, in (c) the Ba surface charac- 
ter and in (d) the tunneling matrix element |M(k, n)| 2 . 
Panel (c) reflects the surface reconstruction, in that the 
Ba projection clearly breaks C4 symmetry. The essen- 
tial point here is that panels (c) and (d) strongly re- 
semble each other, meaning regions of the Fermi surface 



3 



with the largest surface-Ba weight are directly correlated 
with the largest tunneling matrix element. In contrast, 
the weight of the Fe atoms is isotropic around the Fermi 
surfaces (Fig. 2[b]) and does not correlate with the sup- 
pressed tunneling matrix element along the F-Y direction 
evident in panel (d). At low bias voltage, the primary 
tunneling pathway is in fact one from the uppermost Fe 
layer of the terminating As-Fe2-As block to the tip, and 
this runs through the Ba surface atoms, involving the T- 
X region of the Brillouin zone. It is through this filter 
that STS-experiments see the superconducting states of 
these materials. 

Having established the clear correlation between the 
tunneling matrix element and Ba character of the bands, 
we now turn to the calculated tunneling conductance 
dl/dV itself. In Fig. 3 the calculated dl/dV spectra 
obtained for various tip locations above the 2x1 termi- 
nated A- 122 (A=Ba) surface are shown. In simple models 
dl/dV is expected to be a measure of the local DOS of 
the sample at the tip position and so the partial DOS, 
projected on the surface Ba atoms and two representa- 
tive sub-surface Fe atoms is plotted in the upper panel. 
Inspecting the lower panel of the figure, the first observa- 
tion is that the calculated dl/dV spectra exhibit a large 
bias asymmetry, with tunneling into unoccupied states 
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FIG. 2: (Color online) (a) The band structure of the 2x1 
surface reconstructed Ba-122 slab. The Ba weight of the band 
is coded by the thickness of the orange (light) shading. The 
weight of the blue (dark) shading indicates the tip-sample 
tunneling matrix element, (b) the character of the closest-to- 
surface Fe and (c) surface Ba orbitals, respectively, projected 
onto the Fermi surface, (d) A similar Fermi surface projec- 
tion of the tip-sample tunneling matrix element illustrating a 
strong correlation with the Ba projection shown in (c). 
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FIG. 3: (Color online) The tunneling conductance dl/dV 
for the 2x1 reconstructed surface slab, calculated using Eq. 
(1), for various tip positions relative to the crystal lattice. 
The coordinate scheme is sketched in the top right of the 
figure and x = 0, y = corresponds to a tip located directly 
above a Ba atom in the 2x1 termination layer. Upper panel: 
projection of the density of states onto the 3d-orbitals from 
two representative sub-surface Fe (blue, yellow) atoms and 
the 6s/5d-orbitals from the Ba (red) atom of the surface layer. 
Here all dl/dV spectra have been normalized over the energy 
range [0,300] meV. 

of the sample clearly being favored, regardless of the ex- 
act tip position. Thus, consideration of dl/dV rather 
than LDOS now brings the theory into line with the ex- 
perimental data (Fig. 1) in terms of tunneling asymme- 
try. Furthermore, the important intermediary role of the 
adatom states leads to the observation that the main fea- 
tures of the spectra for both electron removal and addi- 
tion correlate well with the structure in the Ba projected 
partial DOS. If it were not for the tunneling matrix el- 
ement in Eq. (1) the simulated dl/dV curve would re- 
semble the LDOS of the Fe atoms, showing - with respect 
to the experimental data - the 'wrong' asymmetry with 
the conductance larger on the occupied, negative sample 
bias side. 

Looking at the three different traces in the lower panel 
of the figure, one sees that the the main features of the 
tunneling spectrum are insensitive to the lateral position 
of the tip position with respect to these high symmetry 
positions within the 2x1 adatom lattice. Experimentally, 
the measured STS spectra are indeed uncorrelated with 
the topographic details within the 2x1 and y/2 x y/2 sur- 
face reconstructions [7] , which lends credibility to the rel- 
evance of both the model and the simulation parameters 
used here. 

In the simulations, for experimentally relevant 
tip^sample separations, the filtering of the tunneling 
current by the adatoms, especially the asymmetry of the 
dl/dV spectra, is also robust for V^xa/2 surfaces. How- 
ever, the peak structure in the calculated dl/dV curves 
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FIG. 4: (Colour online) (a)-(c) DFT-calculated dl/dV maps 
for 2x1 Ba-122 for the energies indicated. A clear contrast- 
switch is observed between 80 and 295 meV positive bias, in 
which the high conductance shifts from being on to between 
the Ba. (d)-(f) Experimental dl/dV maps for a 20 x 20 A 
field of view from a 15% Co-doped Ba-122 crystal (x — 0.3, 
T c = 12 K) at the positive bias voltages indicated. The black 
(white) dashed lines are guides indicating the (missing) Ba 
2x1 rows, (g) A topographic image covering 200 x 200 A. 
The orange square indicates the FOV of panels (d)-(f), the 
topographic data for which is shown in panel (h). 



on the occupied side depends on the assumed locations of 
the sub-surface Fe and As atoms: their positions deter- 
mine the fine details of the surface-related band splittings 
(see online supplementary material). 

In the computer, variation of the tip-sample separa- 
tion goes without risk of a tip crash, and our simulations 
show that dependence of the conductance curves on the 
lateral tip position increases if the tip is moved closer to 
the surface. For the data shown in Fig. 3, a vertical 
tipf^Ba adatom distance of 14A was used, motivated by 
the experimental situation. At tip sample distances well 
below those used typically in experiment (i.e. with the 
tip pushed to within 3 A of the surface), our simulations 
show that the filtering effects through the Ba adatoms 
are reduced. 

Going further than the three different lateral tip posi- 
tions in Fig. 3, we now simulate tunneling conductance 
maps, constructed by varying the tip atom's lateral po- 
sition at a fixed distance from the surface of the slab (14 



A, as before). This generates a set of dl/dV spectra as 
a function of tip position, which in this case are normal- 
ized such that the integral from [0, 300] meV is fixed. To 
create the conductance map for a particular bias voltage 
V, a slice through the simulated data-cube at this bias is 
then plotted as a function of r. 

Applying this procedure to the simulated data for the 
2 x 1-reconstructed surface produces the maps shown in 
panels (a) - (c) of Fig. 4, for V = 80, 150, and 295 meV, 
respectively. At V = 80 meV, the 2x1 reconstruction 
is clearly seen in the calculated maps, appearing - for 
example - as a dark stripe joining the locations of the 
missing Ba atoms. Interestingly, as the bias voltage is 
increased to 150 meV (panel [b]) a clear contrast rever- 
sal takes place. This behavior in the model is confirmed 
in our measured tunneling conductance maps, shown in 
panels (d)-(f) of Fig. 4. These were measured with a 
5 meV modulation on the bias voltage and had set-up 
parameters V set = —300 meV and I set = 50 pA, with 
normalization of the dl/dV spectra to the area of the 
electron addition spectral weight up to 300 meV. The 
(black) white dashed lines in Fig. 4(d)-(f) are guides, 
marking the Ba adatom peaks (valleys) observed in the 
topography (panel [g]), of which a zoom is presented in 
Fig. 4(h). The measured dl/dV maps show a similar rel- 
ative contrast reversal at the same positive bias values. 
The peaks (valleys) in the conductance map at 80 meV 
(panel [d]) reverse roles at 150 meV, appearing instead 
as valleys (peaks) (panel [e]). If one traces along the oc- 
cupied Ba row under the white dashed line in panel [e] 
at 150meV bias, the experiment indicates higher dl/dV 
between the positions of the Ba atoms (the latter visi- 
ble in the topographic scan shown in panel [h] as bright 
features) and the corresponding theory map (panel [b]) 
shows the same behaviour. Moving up to a bias around 
290meV in panel (c), a second type of contrast change 
occurs along the occupied Ba row where the occupied 
(unoccupied) Ba positions become brigher (darker). This 
change is also reflected in the experimental data (panel 
[f]). However, at 295 meV the agreement between theory 
and experiment is not perfect: theory predicts the high- 
est intensity to occur in between the Ba rows whole in 
the experimental data the highest intensity occurs on the 
Ba rows. Nevertheless, the experimental confirmation of 
the contrast changes at lower bias voltage predicted by 
our DFT calculations closes the loop as regards the valid- 
ity and precision of our theoretical description showing a 
heavily Ba-filtered tunneling current in STS of Co-doped 
Ba-122 superconductors. 

In summary, our joint theoretical and experimental in- 
vestigation has uncovered that the alkaline earth adatom 
termination layer in STS experiments on A- 122 iron pnic- 
tide superconductors acts as a tunneling filter. Bands 
with strong Ba5d character in the Ba-122 system - which 
originate in the 2x1 reconstruction - are the door-keeper 
to a sub-set of the Fe3d related states from the upper- 
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most Fe layer of the sub-surface As-Fe2-As block. Our 
model naturally explains the observed strong particle- 
hole asymmetry observed in the STS spectra, which can- 
not be found in the underlying Fe 3d partial DOS. We 
expect very similar adatom-filtering effects to impact the 
tunneling current in other alkaline-earth A- 122 systems, 
but given the sensitivity to the surface reconstruction, 
experimental data regarding the real surface structure - 
such as those from LEED - are required before the dif- 
ferent aspects can be successfully unravelled as we have 
done here. 
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